Astronomy & Astrophysics manuscript no. 1066 


February 2, 2008 


(DOI: will be inserted by hand later) 





Polarization properties of broad absorption line QSOs : new 

statistical clues 

H. Lamy^ and D. Hutsemekers^' 
^ lASB-BIRA, Avenue Circulaire 3, 11 80 Bruxelles, Belgium 

^ Institut d'Astrophysique, Universite de Liege, Sart-Tilman, B-4000 Liege, Belgium 
Received ; accepted 

Abstract. We report the results of several statistical tests performed on a large sample of 139 broad absorption line (BAL) 
QSOs with good quality optical spectra and/or optical polarization data. Correlations between ten optical indices and the 
polarization degree po are systematically searched for. We find six significant non-trivial correlations. In order to identify the 
most important correlations, we perform a principal component analysis with a sample of 30 BAL QSOs and eight quantities 
(including po). Most of the variance (~ 57 %) in the data is contained in two principal components called PCI and PC2. PCI 
is mainly dominated by the correlation between the balnicity index BI and the strength of the Fe II emission; it may be driven 
by the accretion rate of matter onto the central compact object. The variance in PC2 is essentially due to the anti-correlation 
between po and the detachment index DI, indicating that BAL QSOs with P Cygni profiles (DI -C) are usually more polarized 
than those objects with C IV absorption troughs well detached from the corresponding emission lines (DI ^). We show that 
PC2 may be related to the orientation of the BAL QSOs with respect to the line of sight. 

We also present new spectropolarimetric observations of six BAL QSOs. By adding spectropolarimetric data from the literature, 
we build a sample of 21 BAL QSOs for which we define four spectropolarimetric indices describing the polarization properties 
of the absorption and emission lines. We find that the polarization of the C III] emission line is systematically higher than 
the polarization of the C IV emission line, and that the highest polarization in the troughs is correlated to the balnicity index. 
Another important result emerging from the statistical tests performed on this spectropolarimetric sample is a possible anti- 
correlation between the detachment index and a quantity SI which measures the ratio of the depths of the C IV absorption in 
the polarized flux and in the total flux. This correlation indicates that in BAL QSOs with P cygni profiles, the BAL troughs in 
the polarized flux are nearly as deep as in the total flux while, in BAL QSOs with detached absorptions, the BAL troughs in the 
polarized flux are much weaker than in the total flux. 

We show that our main results may be explained in the framework of a "two-component" wind model which is a natural exten- 
sion of the classical wind-from-disk models. In this model, the broad absorption occurs in a dense equatorial wind emerging 
from the accretion disk, while scattering and polarization mainly take place in a polar region. The orientation relative to the 
observer drives the correlations po - DI and DI - SI. While most of our observations can be explained within this framework, 
there are also several indications that other polarization mechanisms, and more particularly resonance scattering, may also be 
at work. 
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1. Introduction 

About 10 to 20% of optically-selected QSOs exhibit broad 
absorption lines (BAL) in the resonance lines of highly ion- 
ized species such as NvA1240, C IV A 1549 or Si IV A 1397 
(Hewett & Foltz EUn^ Tolea et al. lM?^ Reichard et al. 2003 1. 
A smaller fraction 1 %) has in addition broad absorp- 
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tion lines from lower ionization species like Mg II A 2798 or 
Aim A 1858. These objects are often referred to as LoBAL 
or LIBAL QSOs in the literature, in contrast with those ob- 
jects with only high ionization BALs, hereafter denoted HIBAL 
QSOs'. The broad absorption lines are blueshifted relative to 
the emission lines, indicating the presence of outflows with ve- 
locities ranging up to ~ 0.2 c. 

Despite of more than three decades of intense observa- 
tion and modelling, the fundamental nature of the BAL QSOs 

' In order to be consistent with the notations introduced in our pre- 
vious papers, we use here the acronyms HIBAL-LIBAL QSO instead 
of the more widely used HiBAL-LoBAL QSO dichotomy. The differ- 
ences are explained in Sect. 2.2 
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remains controversial. Indeed, the physical mechanism that 
drives the flow is still a matter of debate although there are 
some indications that the outflows are driven by radiation pres- 
sure on spectral lines and originate from an accretion disk ro- 
tating around a massive black hole (Murray et al. 1995 Arav 
119971 Proga et al. 2000 1. Acceleration of the outflow by other 
mechanisms such as gas pressure (Begelman et al. .l991,) . mag- 
netic fields (de Kool & Begelman [1995} or radiation pressure 
by dust (Voit et al. 1993 Scoville & Norman 1995 Everett et 
al. .2002^ has also been proposed and may contribute as well. 
Moreover, it is not yet understood whether BAL QSOs and 
non-BAL QSOs are completely different types of quasars or 
if their properties mainly depend on their orientation towards 
us. 

In order to determine if BAL QSOs and non-BAL QSOs are 
drawn from the same parent population of quasars, Weymann 
et al. ( 119911 hereafter W91) carried out a detailed compara- 
tive study of the emission lines and continuum properties of a 
sample of 42 BAL QSOs and 29 non-BAL QSOs. Recently, 
Reichard et al. (2003 1 made a similar analysis with a much 
larger sample (224 BAL QSOs and 892 non-BAL QSOs) drawn 
from the SDSS EDR (Sloan Digital Sky Survey Early Data 
Release, Stoughton et al. 20021, which essentially confirms 
the results of W91 : there are no significant differences be- 
tween emission line properties of BAL QSOs and non-BAL 
QSOs, except a noticeable enhancement of N V emission in 
BAL QSOs. Also, the continua of the BAL QSOs are redder 
(Reichard et al. I2003> . especially for those objects with low 
ionization absorption lines. This is possibly due to a larger dust 
content (Sprayberry & Foltz ll992l Yamamoto & Vansevicius 
I1999j . Both statistical studies lead to conclusion that BAL 
QSOs and non-BAL QSOs come from a single population. The 
small differences observed in their emission line/continuum 
properties may indicate that BAL QSOs are normal QSOs with 
extreme properties such as larger accretion rates (Reichard et 
al.|2003 Yuan & Wins l2003l . 

Partly based on these results, the most popular paradigm 
considers that a broad absorption line region (BALR) exists in 
all QSOs and has a covering factor ^ 0.1, i.e. close to the rate 
of detection of BAL QSOs among optically-selected QSOs. In 
this orientation scheme, broad absorption lines are observed 
only when our line of sight intercepts the outflow. A com- 
pletely different possibility is to consider that the BAL phe- 
nomenon is an evolutionary phase with significant mass loss 
during the life of a QSO (Voit, Weymann & KoristalW). The 
best argument supporting this evolutionary model is given by 
the large percentage of LIBAL QSOs detected among infrared 
surveys (Low et al. 1989), suggesting that these objects are en- 
shrouded in a cocoon of dust which heavily attenuates the op- 
tical/UV continuum (Sprayberry & Foltz ll99'2t . Moreover, in 
order to explain the very weak [O III] intensity measured in 
those LIBAL QSOs (Boroson & Green 1992 Yuan & Wills 
|2003}, the covering factor of the BALR must be much larger 
than 0.1. 

Since we want to test for orientation as a possible key- 
parameter to understand the BAL QSO/non-BAL QSO di- 
chotomy, polarization may be a useful tool because it is sen- 
sitive to the geometry and orientation of these unresolved ob- 



jects. Previous optical polarization surveys have shown that, as 
a class, BAL QSOs are more polarized than non-BAL QSOs 
and that the polarization is most probably due to scattering 
(Moore & Stockman[T984| Hutsemekers et al l 19981 (hereafter 
Paper I), Schmidt & Hines [T999l Hutsemekers & Lamv l2002l 
Lamy 2003 1. The fact that BAL QSOs are polarized reveals 
that the scattering regions are not spherically symmetric. Many 
current theoretical models of BAL QSOs assume an axisym- 
metric geometry, with either an accretion disk (e.g. Murray et 
al. J_995jt or an opaque dusty torus (e.g. Schmidt & Hines [T997t 
located in the equatorial plane. In this equatorial geometry, the 
BALR emerges either from the accretion disk in the form of a 
continuous flow or it is ablated from the top of the dusty molec- 
ular torus in the form of dense cloudlets. It is then accelerated 
to large velocities. However, there exists recent models which 
locate the BALR within a weak jet in the polar regions in or- 
der to explain the radio-quiet/radio-loud dichotomy and to si- 
multaneously relate the BAL phenomenom to the associated 
absorptions often observed in radio-loud QSOs (Punsly 1 9991 
Kuncic[l999j. In this polar geometry, the BALR is embedded 
in a poorly collimated radio jet and is accelerated by the pres- 
sure of the suiTounding medium which also acts to confine the 
BAL clouds. 

Spectropolarimetry has been used to put constrains on the 
structure of BAL Q SOs (G lenn et al. [T99?l Cohen et al. [T995l 
Good rich & Miller [T9951 Schmidt & Hines [T999l Ogle et al. 
1999). It has revealed the presence of a polarized component 
filling in the absorption troughs and attributed to a continuum 
component polarized by scattering off electrons or dust. The 
scattered and direct rays suffer different amounts of absorption 
which results in a strong increase of the polarization in the BAL 
troughs. Spectropolarimetry has also been used to determine 
the relative size of the different emission/absorption/scattering 
regions in BAL QSOS. However, several axisymetric configu- 
rations are possible to interpret the spectropolarimetric obser- 
vations (Ogle 1997 ). Therefore, we have to find other parame- 
ters which, in connection to the polarization data, may help to 
discriminate among the different models of BAL QSOs. 

In this perspective, we performed in Paper I a systematic 
search of correlations between optical broad-band polarization 
and several indices characterizing the optical spectra of BAL 
QSOs in order to have a clue on the nature of the outflows. 
We find that the main coiTelations involving po could be natu- 
rally interpreted within the "wind-from-disk" (hereafter WfD) 
model of Murray et al. (1995'). However, the size of the sam- 
ple was small (29 objects) and some important but statisti- 
cally marginal trends noticed in Paper I need to be confirmed. 
Schmidt & Hines (1999") also report on a marginal correlation 
between po and the strength of the C IV broad absorption lines, 
which was not detected in Paper I. Since then, a lot of good 
quality polarimetric data of BAL QSOs have been obtained 
(Lamy & Hutsemekers 2000a Sluse et al. 2004 Hutsemekers 
et al. I2004t and the existence of a population of radio loud 
QSOs with intrinsic, BAL-like outflows has been confirmed 
(Becker et al. 2000, 2001 Menou et al.,2^1 ). It is therefore 
important to perform these statistical tests again using a much 
larger sample of BAL QSOs than in Paper I (139 objects in- 
stead of 29), in order to better identify which spectral char- 
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acteristics of BAL QSOs are orientation dependent. This is 
the scope of this paper. The fact that sometimes many prop- 
erties are correlated makes the understanding of the physical 
mechanisms driving them difficult to discern. The Principal 
Component Analysis (PCA) is a statistical method which has 
proven to be very efficient at isolating the most important corre- 
lations in active galactic nuclei (AGN) (e.g. Boroson & Green 
[19921 . Here the PCA will be used with a BAL QSO sample. 

Spectropolarimetry of BAL QSOs with a good signal to 
noise has been published for a large number of BAL QSOs 
(Ogle et al. 1999). In combination with our own spectropolari- 
metric data, we define sp ectropolarimetric indices to perform 
a statistical analysis. The results are then used in combination 
with previous analyses to better constrain the BALR models. 

In Sect. 2, we describe the sets of data. The statistical anal- 
yses are done in Sect. 3. Sect. 4 is devoted to a principal com- 
ponent analysis of the data in order to identify the variables 
that correlate together and to relate them to more fundamen- 
tal underlying parameters that drive the correlations. In Sect. 
5, we present the new spectropolarimetry of 6 BAL QSOs. 
Considering additional data from the literature, spectropolari- 
metric indices are defined and correlations with broad band op- 
tical polarization and with optical indices are searched for. In 
the last section, we discuss and interpret the main results in the 
framework of existing BAL QSO models. 

2. The BAL QSO data 

The sample considered here is made up of 139 BAL QSOs. It is 
essentially a compilation of 93 BAL QSOs coming from large 
polarization surveys. In addition, we also consider 46 BAL 
QSOs with good quality optical (rest-frame UV) spectra but no 
polarization data. The optical polarimetric observations were 
conducted by us (Paper I, Lamy & Hutsemekers 2000a Sluse et 
al. l2004l Hutsemekers et al. i2004) and by Ogle et al. ( 1999 1 and 
Schmidt & Hines ( 1999). The spectra considered for the mea- 
surement of optical indices are those from the W91 and Korista 
et al. J199 3^ samples (mainly drawn from the LBQS), and from 
the recent FIRST (Becker et al. 2000 2001 1 and SDSS/FIRST 
(Menou et al. 2001 1 surveys. Apart from these major surveys, 
we do not consider individually studied BAL QSOs. Indeed 
such objects are not very numerous and are usually observed 
or detected because of their peculiar properties. For the same 
reason, we do not consider the unusual BAL QSOs uncovered 
by Hall et al. (2002). 

The characteristics of the BAL QSOs considered in this 
study are reported in Table 1 and described in detail below. 

2. 1. The broad band polarization data 

The broad band polarization data are compiled from our previ- 
ous surveys (Paper I, Lamy & Hutsemekers I200()al l. from the 
survey of Schmidt & Hines ( 1999^ and from the spectropolari- 
metric atlas of Ogle et al. ( 1999J . When an object is observed 
several times, we choose the value of the polarization degree 
with the lowest uncertainty, i.e. the smallest dp. Our surveys 
have been conducted in the V filter with a typical uncertainty 



~ 0.2 %. These uncertainties are adequate for BAL QSO stud- 
ies since the polarization of these objects peaks around ^ 1 % 
(Paper I, Lamy 2003 1. The measurements of Schmidt & Hines 
(1999), mostly obtained in white light, have often larger un- 
certainties (fjp ~ 0.4%). The data from Ogle et al. (IT9991 
are high-quality broad-band polarizations (4000-8600 A) mea- 
sured from their spectropolarimetric data, and obtained either 
with the Keck or the Palomar telescopes. The typical uncer- 
tainties range from ^ 0.05% for the Keck data to ~ 0.15% 
for the Palomar data. Within the uncertainties, no polarization 
variability is observed for BAL QSOs with multiple observa- 
tions. Finally, a few additional polarization measurements of 
BAL QSOs from the LBQS and FIRST surveys (Becker et al. 
2000 Menou et al.pOOl I have recently been secured using the 
ESO 3.6m telescope at La Silla (Sluse et al. 120041 Hutsemekers 
et al. .2004 ) and are added to the sample. 

The observed degree of polarization p and its uncertainty 
cTp are given in Table 1 together with the debiased degree of 
polarization po according to the Wardle & Kronberg (19741 
method. Since most BAL QSOs in Table 1 have a polarization 
degree of the order of 1 %, we adopt a quality requirement on 
Po for all the statistical tests performed throughout this paper 
and involving this quantity. Data with poor signal-to-noise ratio 
are excluded by considering only those objects with po /<^p > 2 
if Po > 0.80 % or with cTp < 0.40 % if po < 0.80 %. 

2.2. Tlie BAL QSO sub-types 

We use the classification described in detail in Paper I. After a 
careful visual inspection of the optical spectra available in the 
literature, each BAL QSO has been given a number between 2 
and 6, according to the detection and intensity of the low ion- 
ization troughs. HIBAL QSO (type 2) are those objects which 
show broad absorption lines from highly ionized species but do 
not have broad absorption lines from Mg II nor Al III. If Al III 
is not observed in absorption and if the optical spectrum does 
not cover the Mg II spectral region, the BAL QSO is consid- 
ered as an unclassified BAL QSO (type 6). All the objects with 
either Mg II or Al III broad absorption lines in their spectra are 
classified as LIBAL QSOs. According to the strength of the 
low ionization absorption lines, we define three sub-categories 
of LIBAL QSOs : S (strong, type 3), W (weak, type 4) and M 
(marginal, type 5). Note that the term "LoBAL" widely spread 
in the literature essentially corresponds to our type 3. 

2.3. Indices derived from the optical spectra 

We give in Table 1 the following indices originally defined by 
W91 to characterize the absorption/emission features in BAL 
QSOs : (1) the balnicity index BI which is a modified velocity 
equivalent width of the C IV BAL, (2) the half-width at half- 
maximum HWHM and the rest equivalent width REW of the 
C III] A 1909 emission line and of the red part of the C IV emis- 
sion line, (3) the rest equivalent widths of emission features re- 
lated to Fe II and located near 2070 and 2400 A and (4) the 
detachment index DI which measures the onset velocity of the 
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strongest^ C IV BAL trough in units of the adjacent emission 
line half-width, that is the degree of detachment of the absorp- 
tion line relative to the emission one. 

Apart from the W91 sample, none of these indices are given 
in the literature except BI. Therefore, as akeady done in Paper 
I and in Hutsemekers & Lamy (2000 hereafter paper II), we 
have computed these indices whenever possible by using good 
quality spectra digitally scanned from the literature. We fol- 
low the prescriptions given in W9 1 and refer the reader to this 
paper for additional details. The new measurements make use 
of spectra published by W91, Korista et al. J1993I I. Ogle et al. 
( fT9^ . Becker et al. (|2000 2001 1, Menou et al. (|2001 ), and a 
few spectra from the LBQS (Foltz et al. [T989I Chaffee et al. 
[1991] Hewett et al. IT9911 . For one object (B2240-3702), the 
indices have been measured from our spectra (cf. Fig. O. All 
these new measurements are reported in Table 1 together with 
values previously published in Paper I and II. 

In Table 1, following W91, we consider as BAL QSOs 
those objects with a non-zero BI in either C IV or Mg II (since 
usually BI (C IV) > BI (Mg ll)). But for the consistency of the 
statistical analysis, only the BI measured from C IV are consid- 
ered. The BI values in Table 1 are mainly found in the literature 
(W9 1, Kori sta et al. 199 3 Bro therton et al. 1998 2002 Ogle 
et al. [T999l Becker et al. l2000i r2001 Menou et al. 2001 1. For 
the objects considered at the same time in the W91 and the 
Korista samples, we adopt the average of their values. We have 
calculated BI only for those objects observed at the Palomar 
Observatory by Ogle et al. ( 1999> . 

In addition to these indices, we have also reported in Table 
1 the maximum velocity in the C IV BAL trough, Wmax, which 
provides an estimate of the terminal velocity of the flow. For 
the BAL QSOs not included in the samples of Paper I and II, 
the values of Wmax are either given by Becker et al. ( 2000) and 
Menou et al. (2001V or they are evaluated from the Korista et 
al. (19931, Ogle et al. (1999l and Becker et al. (2001 1 spec- 
tra, by measuring, from the blue to the red, the wavelength at 
which the absorption first drops 10% below the flux level de- 
fined by the local continuum. At velocities higher than ^25000 
km s^^, C IV BALs may be contaminated by the Si IV emis- 
sion line, such that measurements of Wmax become inaccurate. 
In a conservative way, we therefore limit t^max to 25000 km 
s^^ from the C IV emission centroid, even if larger values are 
sometimes given by some authors (Becker et al. 120001 Menou 
et al. l2001t . In these cases, 25000 km s~^ constitutes a lower 
limit to the true Umax- Finally, we have measured the slope of 
the continuum, as, by fitting the continuum blueward of C III] 
and assuming a power-law continuum oc i/^"'^ (see Paper 
I for details). 

The errors associated with the measurements of these in- 
dices are difficult to estimate. They are essentially related to 
the positioning of the continuum but they also depend on the 
structure of the lines, the nominal redshift, etc . . . Although the 
uncertainty on the continuum is partly taken into account in its 
definition, errors on BI may still be large. We estimate them 
from values published by different authors (for example the 

^ In objects with multiple absorption troughs, W91 define the 
strongest absorption as the one with the largest equivalent width. 



samples of W91 and Korista et al. 1 19931 1 : typical uncertain- 
ties on BI are around a few hundred kms~^. The typical errors 
on Umax have been estimated from two independent series of 
measurements and are typically around 1000 km s^^. The typ- 
ical error on as is ~ 0.3, as evaluated by varying the possible 
continua (Paper I). The error associated with the detachment 
index DI is typically around 15 % from a comparison between 
our values and those from W91, although it could be higher in 
some cases, especially for those objects with multiple troughs 
or complex absorptions. Indeed, depending on the continuum, 
one might select different troughs as being the strongest one. 
By doing the same kind of comparison with W91, we estimate 
the eiTors associated with C IV HWHM and C IV HREW to be 
10% and 20% respectively . 

Finally, we have included in Table 1 the B absolute mag- 
nitude of the objects, Mb, which is calculated by follow- 
ing the prescriptions of Veron & Veron |2000it (with Hq = 
50kms^^ Mpc~^, qo = and an optical spectral index a = 
0.3) and by applying the "BAL K-correction" of Stocke et al. 
(1992). 

3. Analysis of correlations involving the 
polarization and spectral indices 

Possible coiTelations involving the optical polarization and the 
indices defined in Sect. 2.3 have been searched for by comput- 
ing the Kendall r and the Spearman rank correlation coeffi- 
cients (Press et al. .l989j . While our main goal is to investigate 
the behavior of the polarization, correlations between the in- 
dices themselves have also been searched for since several new 
measurements are available. 

The results of these correlation tests are given in Table El 
The first two columns list the two indices we test. The next two 
columns give the probabilities P^- and Pr^ to find a coiTelation 
by chance between two uncorrected quantities. The column la- 
beled n is the number of objects involved in the analysis. Since 
we search for correlations among 1 1 quantities, a total of 55 
correlations tests are performed. Since the number of objects 
is rather large for most pairs of indices, we may reasonably 
consider as significant only those correlations with P < 0.02. 
With 55 correlations tested, we expect to detect one false corre- 
lation with such a probability, assuming they are independent. 
In general, an excellent agreement is found between the results 
of the two statistical tests. The correlations have also been in- 
vestigated for the HIBAL and LIBAL QSO subsamples sepa- 
rately and the results are given in Table |2l A few correlations 
are detected only when considering the whole BAL QSO sam- 
ple, illustrating the importance of the size of the sample. 

Among all significant correlations reported in Table|2j sev- 
eral link in an obvious manner quantities either related by their 
definition or by a trivial common physical origin. In the first 
category, as an example, the correlation between BI and Umax 
simply reflects the fact that the wider the trough the larger 
both BI and Umax (with BI < Umax " 5000 according to their 
definitions). Another such trivial correlation is DI - Umax^ the 
most detached absorptions have necessarily large Umax- In the 
second category, we may quote the coiTelations between C IV 
HWHM and Cm] HWHM and between the Fell indices. 
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Table 2. Analysis of correlations between po and the indices described in Sect. 2.3 as well as correlations between the indices 
themselves. Three samples of BAL QSOs are considered, n is the number of objects considered in the tests. The last column 
indicates whether the correlation is positive or negative (identical for the three BAL QSO samples). Boldfaced entries are the 
least likely to be due to random chance. 







BAL QSOs 




HI BAL QSOs 


LI BAL QSOs 




Index 1 


Index2 






n 


Pr 


Prs 


n 


Pr 


Prs 


n 




Po 


rsl 


A 1 'TiC 

0.176 


A 1 OA 

0.180 


<o 
5o 


0.557 


A £A£ 

0.646 


OT 

27 


A 1 CO 

0.153 


A 1 OA 

0.180 


O 1 

21 


+ 


Po 


DI 


0.001 


A AA1 

0.003 


54 


0.437 


0.469 


oc 

25 


A A A /I 

0.004 


A AAO 

0.003 


20 




Po 


C IV HWHM 


0.834 


0.824 


54 


0.772 


0.890 


26 


0.888 


0.898 


19 




Po 


C IV HREW 


A 1 1 '5 

0.113 


A AAA 

0.099 


eo 

53 


A COO 

0.538 


A C/CA 

0.560 


oc 

25 


A f\10 

0.078 


A nAdZ 

0.046 


1 A 

19 


+ 


Po 


/"I nil 

L. IIIJ HWHM 


A ACA 

0.059 


A AOC 

0.085 


A O 

48 


A nA A 

0.744 


A OCA 

0.859 


O A 

24 


A OOO 

0.322 


A .1 CO 

0.458 


16 




Po 


C IIIJ REW 


A OAC 

0.305 


A OAC 

0.305 


oo 
33 


0.763 


0.770 


1 c 

15 


A 1 nf\ 

0.170 


A OAO 

0.308 


1 o 

12 


+ 


Po 


re II A z4UU 


0.41 1 


0.397 


O 1 


0.227 


0.217 


17 


0.663 


0.792 


13 


+ 


Po 


re II A 2070 


0.022 


0.033 


38 


0.007 


0.003 


17 


0.516 


0.541 


15 


+ 


Po 


as 


A ACI 


A A0 1 


.1 A 

49 


A AOA 


A AO /I 


oo 


ACT? 
U.51 / 


A AC\A 

U.494 


1 o 
18 


+ 


Po 


I'max 


0.978 


0.967 


J 1 


0.656 


0.631 


27 


0.710 


0.627 


21 


+ 


BI 


DI 


0.208 


0.213 


93 


0.022 


0.036 


43 


0.390 


0.285 


26 


+ 


BI 


C IV HWHM 


A T/l 1 

U. /4i 


U. /SU 


92 


U.469 


A /I C 1 

U.45 i 


A O 

43 


A C 1 O 

U.513 


A COO 

U.5o3 


oc 

25 


+ 


rsl 


Civ HREW 


0.495 


0.405 


AA 
9U 


A AiCC 

0.965 


0.869 


42 


0.301 


0.209 


OC 

25 




BI 


L- IIIJ HWHM 


0.044 


A A'5A 

0.030 


iCA 

69 


A 1 AA 

0.100 


A AAO 

0.093 


A 1 

41 


A oo /( 

0.384 


A Ann 

0.409 


1 o 

18 


+ 


rsl 


C IIIJ REW 


A CAC 

U.51J5 


A A 
U.466 


o o 
3o 


A OOO 
(J. 080 


A A 1 O 

U.912 


1 A 

19 


A OTO 

U.2/3 


A O C/1 

U.354 


1 o 

Iz 


+ 


rsl 


TT \ OATA 

re II AzO*0 


A AAA 
0.000 


A AAA 
0.000 


43 


0.103 


0.062 


21 


0.067 


0.087 


15 


+ 


BI 


bell A 2400 


A AAA 
0.000 


A AAA 

0.000 


35 


A AAO 

0.008 


A AA«C 

0.006 


21 


0.197 


0.254 


13 


+ 


rsl 


as 


A AA1 

O.OOj 


A AA1 
O.OOJ 


66 


A nr\Q 


A HA 1 
U. /41 


oo 
32 


A AO 1 


A A1 £ 
O.OlO 


Oyl 
Z4 


+ 


rsl 


'I'max 


A AAA 

0.000 


A AAA 

0.000 


1 AO 

lUo 


A AAA 

0.000 


A AAA 

0.000 


CO 

53 


A AOA 

0.030 


A ACO 

0.058 


OA 

30 


+ 


DI 


C IV HWHM 


0.527 


0.509 


91 


0.200 


0.224 


42 


0.513 


0.619 


25 




Dl 


C IV HREW 


A AA'^ 

0.002 


A AA1 
0.001 


AA 
9U 


A A A /I 

0.004 


A AA1 
0.001 


41 


A AA^ 

0.006 


A fkfiA 

0.004 


25 




DI 


C IIlJ HWHM 


0.247 


0.217 


65 


0.053 


0.047 


38 


0.970 


0.798 


18 


+ 


DI 


ClIlJ REW 


0.040 


0.054 


38 


A CTC 

0.575 


0.459 


19 


0.217 


0.138 


1 o 

12 




Dl 


re II AzO/0 


A OK 

U.335 


A '3 1 A 


A O 

4j 


A OAC 

U.3U5 


A OAT 


1 

zl 


A COA 


U.5 /6 


1 c 

15 


+ 


Dl 


XT \ o Ar\r\ 
re II A 2400 


0.881 


0.092 


o c 
35 


0.398 


0.336 


21 


0.243 


0.343 


13 


+ 


DI 


as 


0.532 


0.469 


61 


0.836 


0.868 


29 


0.843 


0.922 


22 




DI 


I'max 


A AAA 
0.000 


A AAA 
0.000 


AO 

93 


A AAA 
0.000 


A AAA 
0.000 


A O 

43 


A AylA 
U.U4U 


A AOO 


o<c 
26 


+ 


C IV HWHM 


CiV HREW 


0.387 


n At c 

0.415 


AA 
9U 


0.616 


A /no 

0.678 


/(O 

42 


A AOO 

0.082 


A 1 O 1 

0.131 


oc 

25 


+ 


C IV HWHM 


C IIIJ HWHM 


A AA/I 

0.004 


A AAC 

0.005 


iZA 

64 


A 1 AC 

0.105 


f\ t AO 

0.148 


0'7 

37 


A fk-t "T 

0.017 


A A1 1 

0.011 


1 o 

18 


+ 


C IV HWHM 


L. IIIJ REW 


U.54D 


A COA 

U.569 


in 
i 1 


A AAA 
U.9U9 


A AO C 

U.925 


1 o 

lo 


1 AAA 
I .UUU 


A AO 1 

U.93i 


1 o 

iz 




C IV HWHM 


TT \ OATA 

re II A 20/0 


0.204 


0.205 


A 1 

41 


0.192 


0.229 


20 


0.702 


0.887 


1 A 

14 


+ 


C IV HWHM 


bell A 2400 


A iCOA 

0.630 


A /CO C 

0.635 


oo 

33 


A ACA 

0.059 


A AC 1 

0.051 


OA 

20 


A 1 O 1 

0.131 


A 1 AA 

0.199 


1 o 

12 


+ 


U IV HWHM 




A /I /n 
U.44Z 


A Q OA 


C A 

59 


A O CO 

U.o5o 


A O /f O 


oo 

zo 


A AOA 


A AO O 

U.U35 


O 1 

zl 


+ 


L- IV HWHM 


^max 


A A1 A 
0.014 


A A1 C 
0.015 


AO 

9z 


0.066 


0.094 


43 


0.228 


0.221 


25 


+ 


C IV HRBW 


C IIlJ HWHM 


A AA1 
0.003 


A Al 1 
0.011 


63 


A AA^ 

0.006 


A A1 O 

0.012 


36 


0.340 


0.373 


18 




C IV HREW 


CIIIJ REW 


A AAA 

0.000 


A AAA 

0.000 


3 / 


0.027 


0.029 


18 


0.095 


0.042 


12 


+ 


C IV HREW 


H7« XX \ OA'7A 

Fell A 2070 


A Ado 

0.498 


0.525 


A 1 

41 


A 0.1 C 

0.845 


A OOC 

0.885 


OA 

20 


0.868 


A nn£ 

0.976 


1 A 

14 




C IV HREW 


re II A 2400 


0.938 


0.985 


33 


0.896 


0.830 


20 


0.403 


0.555 


1 o 

12 




C IV HREW 




A CAT 


A ^QA 


CO 

5o 


A AOO 

U.953 


A OOO 

U.80Z 


OT 


A OTO 

U.3 /o 


A OOA 

U.339 


O 1 

Zl 




U IV HREW 


^max 


0.069 


0.042 


AA 

9U 


0.029 


0.028 


42 


0.905 


0.830 


25 




C III] HWHM 


C III] REW 


0.512 


0.498 


38 


0.916 


0.932 


19 


0.272 


0.183 


12 


+ 


L IIIJ HWHM 


Fell A 2070 


U.2/5 


A '^^O 
U.ZZS 


A 1 

41 


A ACO 

U.952 


A O^T 

U.o6/ 


1 

zl 


A AOO 


A A^l 

U.U61 


1 o 
13 


+ 


U IIIJ HWHM 


II \ OA A A 

re II A z4UU 


0.852 


0.867 


oo 
33 


0.504 


0.482 


21 


0.392 


0.247 


1 1 


+ 


C III] HWHM 


as 


0.039 


0.046 


56 


0.282 


0.312 


30 


0.084 


0.171 


17 


+ 


C III] HWHM 




050 


06^ 


69 


1 56 


9^7 


41 


090 


099 


1 8 


+ 


C III] REW 


Fp TT \ 9070 
rc 11 A zu 1 u 


0.920 


0.912 


DO 


0.574 


0.557 


19 


0.411 


0.457 


12 




C III] REW 


Fell A 2400 


0.943 


0.997 


30 


0.207 


0.193 


19 


0.025 


0.038 


10 




Cm] REW 


as 


0.880 


0.967 


38 


0.674 


0.881 


19 


0.681 


0.527 


12 


+ 


C III] REW 




0.125 


0.117 


38 


0.379 


0.501 


19 


0.726 


0.678 


12 




Fell A 2070 


Fell A 2400 


0.000 


0.000 


35 


0.040 


0.026 


21 


0.057 


0.079 


13 


+ 


Fell A 2070 


as 


0.133 


0.110 


43 


0.763 


0.771 


21 


0.347 


0.226 


15 


+ 


Fell A 2070 


Vmax 


0.005 


0.003 


43 


0.414 


0.383 


21 


0.096 


0.080 


15 


+ 


Fell A 2400 


as 


0.048 


0.079 


35 


0.116 


0.167 


21 


0.581 


0.648 


13 


+ 


Fen A 2400 


''max 


0.001 


0.000 


35 


0.057 


0.038 


21 


0.570 


0.406 


L3 


+ 


an 


f'max 


0.001 


0.001 


66 


0.280 


0.226 


32 


0.008 


0.005 


24 


+ 



6 



Lamy & Hutsemekers; Polarization properties of broad absorption line QSOs : new statistical clues 



□ m ■ □ X 
ok 



: X n n'^ 

m x„ 

-i □■ □ □ 

1 □ ^ 

- 1 /^xft d*x 



□ 

□ □ ^ 

^ , \ 

20 





mo ^ 


□ 






m ■ 




X ^ 




□ □ 


□ 


■ 


□ ^ 


■ 





0.0 0-5 



Fig. 1. The most significant non-trivial correlations (see text), pa is the debiased broad-band optical polarization (in %), DI is 
the detachment index (unitless), C IV HREW is the rest equivalent width of the red part of the C IV BEL (in A), BI is the 
balnicity index (in 10'^ km s~^), as is the slope of the continuum measured shortward of C III] (unitless) and Fe II 2400 is the 
rest equivalent width of the Fe II emission feature near 2400 A. Open squares are HI BAL QSOs, filled squares are LI BAL 
QSOs and crosses are unclassified BAL QSOs. 



Finally, a few correlations in Table |2l are due to the fact that 
two quantities are correlated with a third one : Wmax is corre- 
lated with Fe II and a b because, as we will see below, the bal- 
nicity index is correlated with these quantities. With all these 
arguments in mind, we are left with six significant non-trivial 
correlations : po - DI, BI - Fe II A 2400, BI - a^, DI - C IV 
HREW, pq - Fe II a 2070 and C IV HWHM - w^ax- They are 
all discussed in detail below and the most interesting ones are 
illustrated in Fig. [2 We have checked that these coiTelations 
remain significant when we exclude the objects with the most 
extreme properties (e.g. the two LIBAL QSOs with DI > 10 in 
Fig.CJi). 

The degree of optical polarization, pq, is essentially uncor- 
related with most indices apart from a few exceptions. The 
most striking one is the anti-correlation between po and the 
detachment index DI, already detected in Paper I and first sug- 
gested by Goodrich J1997> . It is confirmed here with a higher 
significance (P ~ 0.002) and on the basis of a larger sample. 
It means that, in average, BAL QSOs with C IV P Cygni pro- 



file ^ are more polarized than objects with detached absorption 
troughs. As can be inferred from Table |2] this anti-correlation 
is strong for the whole BAL and the LIBAL QSO samples, but 
not detected in the HIBAL QSO sample. Fig.^shows that the 
relation between p^ and DI looks more like an inequality than 
a real anti-coiTelation. Indeed, for small values of DI, pq varies 
within a large range of values (from to 5 %), while for large 
values of DI, po is smaller. This coiTelation indicates that DI 
may play a key-role in understanding the BAL QSO physics, as 
suspected by Hartig & Baldwin ( 1986 1. Since it is the most im- 
portant correlation involving polarization, it will be discussed 
in detail in Sect. 6. 

The correlation between po and Fe II A 2070 is more doubt- 
ful because it is significant only for the HIBAL QSO subsample 
(marginal for the whole BAL QSO sample). Moreover, the fact 
that Po increases with the iron index does not support the usual 
interpretation of a dilution of po by an unpolarized pseudo- 



^ DI is small in P Cygni type profiles since the absorption starts near 
zero velocity 
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continuum formed of a large number of metastable iron emis- 
sion lines. Finally, we do not see any obvious reason why po 
should correlate with Fe II A 2070 and not with Fe II A 2400. 
For all these reasons, we believe that this correlation is proba- 
bly due to chance. 

It is also important to remark that we do not detect any cor- 
relation between po and BI. Such a correlation was reported by 
Schmidt & Hines ( 1999 1 in their Fig. 8 (although their corre- 
lation was marginal with P ^ 0.05). Neither do we recover 
the correlation between po and the slope of the continuum a b 
discussed in Paper I and already interpreted as an artifact. 

The correlations between BI and the Fe II A 2400 index on 
one hand and between DI and C IV HREW on the other hand 
have been reported previously in W9 1 . They are confirmed here 
with larger samples and with the use of BI averaged from the 
measurements of W91 and Korista et al. U993i) . The corre- 
lation between BI and the Fe II A 2400 index remains strong 
within the HIBAL QSO subsample while it disappears in the 
LIBAL QSOs subsample. These two correlations are illustrated 
in Fig. n Another significant correlation involving BI is the 
correlation between BI and the slope of the continuum as (see 
Fig [D- It was found in Paper I and is confirmed here. It indi- 
cates that BAL QSOs with redder spectra display on average 
stronger absorption lines. Since these two correlations involv- 
ing BI appear to be important in the description of the BAL 
phenomenon, they will be discussed later in the framework of 
the BALR models (Sect. 6.4). 

The anti-correlation between DI and C IV HREW may 
possibly be explained in terms of radiative pumping of the 
BEL photons. Indeed, for P Cygni profiles (DI <), the BAL 
trough cuts off the blue wing of the C IV broad emission line 
(Turnshek ll988> . Photons are absorbed and partly re-emitted in 
the red wing, therefore increasing the equivalent width of the 
red part of the C IV emission. 

Finally, the correlation between C IV HWHM and Umax 
—definitely less significant— is not obvious to interpret. 
Although it may be spurious, it may also suggest that the dy- 
namics of the BAL wind is connected in some way to the dy- 
namics of the broad emission line region from which it could 
be launched. 

4. Principal Component Analysis 

When we are dealing with a large number n of BAL QSOs and 
a large number p of distinct parameters characterizing the QSO 
properties, the study of the correlations two by two may be 
quite complex. Indeed, as can be seen from Table|2] several in- 
dices appear correlated, which makes it difficult to disentangle 
the leading correlations and consequently the physical mecha- 
nisms behind them. The Principal Component Analysis (PCA) 
is an alternative mathematical tool specifically designed to deal 
with such multivariate problems. Basically, the PCA identifies 
those parameters that correlate together and decreases the com- 
plexity of the problem by gathering them into only a few new 
parameters. In a p-dimensional space, PCA identifies those di- 
rections in which the cloud of points representing the QSOs is 
most elongated and uses them to define a new system of coordi- 
nates. Each axis in this new system is called an eigenvector and 



is simply a linear combination of the initial parameters. The 
first eigenvectors then contains most of the variance observed 
in the sample and identify the most significant correlations. For 
a complete mathematical description of multivariate statistics, 
see Murtagh & Heck J1987> . For nice illustrations of the ap- 
plication of PCA to some AGN data samples, see Boroson & 
Green ( il992ii . Francis & Wills ( fT999b and Boroson (,2002< . 

Table 3. Results of the Principal Component Analysis. The first 
three principal components out of a total of 8 are given in order 
of their contribution to the total variance. The first row gives 
the variances of the data along the direction of the correspond- 
ing principal component. The sum of all the variances is equal 
to 8, the number of variables considered in the analysis. The 
variance of each principal component is given as "proportion" 
in the second row and as "the cumulative proportion" in the 
third one. The rest of the table gives the correlation coefficients 
between the original variables and the principal components. 
Significant correlations are boldfaced, i.e. those ones ith coef- 
ficients whose absolute values are > 1/ \/8 = 0.3536 
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0.4935 


-0.1053 


0.3406 



Since some objects have different wavelength coverages 
and/or their spectra are of poor quality, there are obviously 
missing data for many QSOs of our sample. Because a full 
data set is needed to run the PCA, the number of useful ob- 
jects will then depend on the indices taken into account. As 
a rule of thumb, the more numerous the indices, the smaller 
the sample. Moreover, we wish to avoid most of the obvious 
correlations between closely related variables which would be 
dominant in the first eigenvectors. After several trials, we made 
a compromise between the size of the sample and the number 
of the adopted parameters. Together with pq, all the parame- 
ters describing the emission and the absorption properties of 
the CiV line profile (BL DL C IV HREW, C IV HWHM and 
^^max) are considered. Indeed, C IV measurements are available 
from medium resolution spectra for a large number of QSOs of 
our sample (W91, Korista et al. I1993I I. We also consider one 
of the Fe II intensity indices : Fe II A2400. The Fe II A2400 in- 
dex is favored since it is much stronger than the A2070 index. 
Finally, we add the slope of the continuum, as, another im- 
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portant quantity in the description of BAL QSO spectral prop- 
erties. The indices related to C III] are not considered because 
of the well-known blend with Al III and a possible blend with 
Fe III in some LIBAL QSOs (Hartig & Baldwin 1986 1. Several 
PCA tests were performed including the absolute B magnitude 
of the objects. Since this quantity does not correlate with any 
of the indices, it may be safely discarded from the analysis. 
However, we cannot definitely rule out any correlation with 
the intrinsic luminosity since some BAL QSOs may be signif- 
icantly reddened due to extinction (Sprayberry & Foltz 1992| 
Reichard et al. l2003> . The final sample then consists in 30 BAL 
QSOs with good quality measurements for the 8 variables pre- 
viously described. 

The PCA was carried out with the Fortran code written by 
Murtagh & Heck (19871. Since the ranges of the variables are 
very different, the original data are first centred and reduced to 
unit standard deviation (correlation matrix option). 

We present the results in Table|3]which lists the three most 
significant principal components (PC) and their projections 
upon the eight variables described above. Only the PCs with 
a variance (eigenvalue) larger than 1 are considered. More than 
half the total variance of the sample (57 %) is accounted for by 
the first two PCs which are elongated with variances larger than 
or of the order of twice that of the original data. We concentrate 
therefore exclusively on these two principal components which 
we call PCI and PC2. 

Apart from the usual redundant correlations, PCI is dom- 
inated by the correlation between the balnicity index BI and 
the strength of the Fe II emission (and marginally as), while 
PC2 correlates essentially po, DI and C IV HREW. This anal- 
ysis is in good agreement with the classical correlation tests 
performed in Sect. 3. Indeed, two of the four important cor- 
relations discussed in Sect. 3 (po — DI and DI — C IV 
HREW) dominate the eigenvector PC2 while a third one (BI 
— Fe II A 2400) dominates PCL Small differences with the pre- 
vious analysis (i.e. the significance of some correlations) may 
be explained in terms of the smaller sample considered here. 
Therefore, the main result emerging from the PCA is that BI, 
Fe II A2400 and maybe as on one hand, and po, DI and C IV 
HREW, on the other hand, vary together. As will be discussed 
in Sect. 6, PCI could be connected to the accretion rate of mat- 
ter (since this interpretation may qualitatively reproduce the 
correlation BI - Fe ll), while PC2 could be related to the ori- 
entation of the BAL QSO with respect to the line of sight. The 
existence of two PCs outlines that, although important, orien- 
tation is probably not the only parameter governing the BAL 
QSO properties. 

Fig. El illustrates the distribution of the 30 BAL QSOs in 
the PC1-PC2 space. Different symbols have been used to sep- 
arate the HIBAL and LIBAL QSO subsamples and to distin- 
guish the lower luminosity objects (with Mb > —28) from the 
brighter ones (with Mb < -28). The LIBAL QSO subsample 
is essentially located in the PCI < region while most of the 
HIBAL QSOs have PCI > 0. The HIBAL and LIBAL QSO 
subsamples are therefore better discriminated with PCI than 
with PC2. Indeed, the LIBAL QSOs have the largest BI and 
Fe II emission or at least they display the most extreme values 
(see Fig.Q. If PC2 is indeed related to orientation, it plays a 
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Fig. 2. Distribution of the 30 BAL QSOs with respect to PCI 
and PC2. Open squares represent HI BAL QSOs, filled squares 
LI BAL QSOs while the cross represents the only unclassi- 
fied BAL QSO considered in the PCA. Large symbols refer 
to bright QSOs with Mb < —28 and small symbols represent 
fainter QSOs with Mb > -28. 

less important role in explaining the differences observed be- 
tween HIBAL and LIBAL QSOs. Among the LIBAL QSOs, 
there is also a trend (not represented in Fig.|2ji that the objects 
with stronger low ionization troughs have more negative values 
of PCI, i.e. the S LIBAL QSOs have more negative PCI than 
W or M LIBAL QSOs, the latter having values of PCI closer 
to those of HIBAL QSOs. Finally, interestingly enough, fainter 
QSOs appear to be located in the PCI < part of the diagram. 

5. BAL QSOs spectropolarimetry 

5. 1. Description of tine observed sample 

Seven BAL QSOs were chosen for spectropolarimetry because 
of their brightness, their high degree of linear polarization in 
the V filter (Paper I) and their redshift adequate for the obser- 
vation of the C IV absorption line. Moreover, all these objects 
are LIBAL QSOs with strong, weak or marginal absorptions 
from low ionization species. They also span a range of other 
properties : B2225-0534 and B2240-3702 are BAL QSOs 
with P Cygni CiV profiles (DI <) while B1246-0542 and 
B 101 1+0906 are objects with detached C IV troughs (DI >). 
B1413+1143 is a gravitationally lensed QSO with evidence 
for a micro-lensing effect (Angonin et al. 119901 Hutsemekers 
[T993t and B0059-2735 is the prototype of the so-called "Iron 
Lo BAL QSOs" (Hazard et al. 1987. Becker et al. |T997}. 
Spectropolarimetry of B0059— 2735 has been pubUshed else- 
where (Lamy & Hutsemekers l2000bt and will not be discussed 
again. 

Spectropolarimetry of B2225— 0534 has been abundantly 
described in the literature in the last two decades (Stockman 
et al. 1981. Cohen et al. .l995 Goodrich & Miller J995. Ogle 
et al. I1999> . Since there is no indication of variabiUty in this 
object, these data have enabled us to check the whole observa- 
tion/reduction procedure. 
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Fig, 3. Spectropolarimetric data for the BAL QSOs B1246-0542 (top left), BlOl 1+0906 (top right), B2225-0534 (bottom left) 
and B2240— 3702 (bottom right). For each object, there are six panels; from top to bottom : total flux F\, degree of polarization 
p, associated error ap, polarization position angle 9, associated error ag , and polarized flux p x Fx.p has not been debiased. The 
main absorption lines are indicated. 
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Fig. 4. Spectropolarimetric data for the BAL QSOs B1331-0108 (left) and B1413+1143 (right). Panels are as in Fig.|3] The 
B300 and R300 spectra were cut on their reddest and bluest parts respectively and merged together at the wavelength indicated 
by the vertical dashed line. The main absorption lines are indicated as well as features due to cosmic rays hits. The asterisk in 
the spectrum of B1413+1143 indicates an intervening absorption system at z ^ 1.66. The Earth symbol indicates a telluric 
absorption. 



5.2. Observations and data reduction 

The spectropolarimetric observations were carried out at the 
3.6m telescope (La Silla, ESO, Chile) equipped with EFOSCl 
during two runs in April 5-8, 1995 and September 10-12, 
1996. 

Spectropolarimetry is performed with EFOSCl by insert- 
ing a Wollaston prism in the filter wheel and a grism in the 
grism wheel (di Serego Alighieri ll998l l. The Wollaston prism 
splits the incoming radiation into two orthogonally polarized 
spectra separated by 20'.' Each object positioned on the slit has 
two spectra projected on the CCD detector. The measurement 
of their ratio at two different Wollaston angle positions sepa- 
rated by 45° gives the two Stokes parameters describing the 
polarization (no half- wave plate was available at that time). 
In practice, the Wollaston prism does not move but rotates to- 
gether with the whole EFOSCl instrument, including the slit. 
In our case, this is not critical since our objects are point-like. 
In order to avoid overlapping of the object and sky spectra, 
a mask is placed in front of the slit; it is transparent only in 
slots whose length is equal to the beam separation produced by 
the Wollaston prism. The 2D image recorded on the CCD then 
consists of a number of alternate bands, two of them containing 



the two orthogonally polarized spectra of the object, the others 
being devoted to the sky. 

The detector was a 512 x 512 Tek CCD (ESO#26) with 
an inverse gain of 3.8 electrons per ADU and a read-out-noise 
of 7.9 electrons. The pixel size was 27 /im corresponding to 
0'.'605 on the sky (Melnick et al. 1989). The sUt had a width 
of 2'.' The seeing was typically around 1'.'4. The B300 grism 
was used for the seven objects giving a spectral coverage of 
3600-6800 A. For two objects, the R300 grism was also used to 
enlarge the useful spectral range up to ^ 8300 A allowing us to 
observe the Mg II BAL. The spectral resolution was 6.7 A per 
pixel. For each grism and position angle of the Wollaston, three 
frames were secured for a total exposure time of 90 minutes per 
Wollaston position. The splitting in shorter exposures ensures 
that the contamination by cosmic rays can be corrected. 

Data reduction was accomplished with procedures devel- 
oped within the ESO MIDAS package. The 2D raw data are 
first corrected for bias and flat-field. Cosmic ray events are 
interactively corrected and replaced by the median value of a 
small surrounding area. However, it is sometimes impossible 
to correct for some cosmic rays falling on the object spec- 
tra. These cosmic rays are easily detected by comparing the 
individual frames and they are flagged in the final ID spec- 
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tra with a 'hit' label. The two object spectra and the four sky 
spectra (two on both sides of the object) are individually ex- 
tracted with exactly the same procedure. The six spectra are 
calibrated in wavelength using He+Ar internal lamp spectra. 
Since we observed small offsets between corresponding He+Ar 
emission lines in orthogonally polarized spectra, six individ- 
ual wavelength calibration curves were produced, and each ob- 
ject/sky spectra was independently wavelength calibrated. The 
sky spectra are averaged on both sides of the object and sub- 
tracted from the object spectra. Spectra are then corrected for 
the atmospheric extinction. Since there was no HWP plate in 
the setup, a flux calibration was done in order to correct for 
a possible difference of transmission of the Wollaston prism + 
grism for the ordinary and extraordinary rays. For that purpose, 
we observed spectrophotometric standard stars, assumed to be 
unpolarized, with the same instrumental setup. The standard 
star EG274 (Stone & Baldwin|T983 Baldwin & Stone I984i) 
was found to be very stable from night to night and was used to 
obtain the individual flux calibration curves. For that star, we 
measure a weak polarization degree (p < 0.3%), more or less 
constant with wavelength, and in agreement with the observa- 
tions by di Serego Alighieri (il994i> with the same instrumen- 
tal setup. Although the instrumental polarization induced by 
the instrument is therefore small, we took account of it by per- 
forming the flux calibration. Finally, it was necessary to rebin 
the data in wavelength in order to achieve a polarization uncer- 
tainty of ^ 0.2 - 0.3 % in the continuum. The spectra are re- 
binned on two original pixels, which corresponds to a spectral 
resolution of about 13 A. The associated errors are calculated 
by propagating the errors from the photon noise in the object 
and sky spectra. 

The Stokes parameters q and u, the degree of polarization p 
and the angle of polarization 9 are computed from these re- 
binned spectra in the usual way (di Serego Alighieri fTWi. 
Rebinning is done before the calculation of the Stokes parame- 
ters since q and u have non-gaussian distribution errors (Clarke 
& Stewart'1986>. The photon noise is propagated throughout 
the equations in order to obtain cTp and ere in every bin of wave- 
length. Although p is a well-known biased quantity (Simmons 
& Stewart i 19858 . we do not de-bias it since our data have large 
S/N ratios. Finally, we calculate the polarized flux by multiply- 
ing the total flux by the degree of polarization. All these results 
are given in Figs. |3l and |3 for the six BAL QSOs . 

A polarimetric standard star, HD 16 1291, was observed dur- 
ing both runs with the same setup. The interstellar polariza- 
tion wavelength dependence was found in good agreement with 
that measured by di Serego Alighieri et al. (J994j using the 
same setup, and with the parameters given by Serkowski et al. 
ltT975t . 

5.3. Description of individual objects 

The spectropolarimetry of individual BAL QSOs is discussed 
in this section. We emphasize mainly their differences as well 
as new results that were not previously reported by Schmidt 
& Hines (1999' I and more particularly by Ogle et al. (1999) 
in their spectropolarimetric atlas. Note that Ogle et al. (il999i) 



use q' instead of p. q' is the Stokes parameter q rotated by the 
angle 2 6 where 9 is the mean continuum polarization angle. 
Although q' is unbiased, it has the disavantadge that its value is 
difficult to interpret when measured in emission or absorption 
lines with significant rotation of 9. 

1. B1246-0542 This BAL QSO has deep broad ab- 
sorption lines clearly detached from the coiTesponding emis- 
sion lines. The polarization rises strongly in the C IV and Si IV 
broad absorption Unes, mainly in the deepest troughs, while it 
slightly decreases in the emission lines. There is no rotation 
of the polarization accross those lines. The polarized flux of 
B 1246— 0542 is remarkable: first, it is nearly featureless with 
no trace of absorption, which is quite unusual among BAL 
QSOs (Schmidt & Hines .1999. Ogle et al. 1999) . Second, it 
shows an unusually large N V emission, with no other emis- 
sion lines. Apparently unnoticed, such a large N V emission is 
in fact observed in the polarized flux of several BAL QSOs 
from the spectropolarimetric atlas of Ogle J1999> . as well as 
in B 101 1+0906. A possible interpretation of this higher Nv 
emission in the polarized flux, in contrast with other lines such 
as C IV and Si IV, could be an extra N V emission due to reso- 
nance scattering of Lya photons (Surdej & Hutsemekers I 1 987) . 
polarized due to a non-spherical geometry of the BALR (Lee 
& Blandford [T997l . 

2. B 101 1+0906 The spectrum of this BAL QSO 

has some similarities with B1246— 0542, with deep absorp- 
tions well detached from relatively faint emission lines, and 
a strong N V emission. As usual, the polarization increases 
within the C IV trough. The polarized flux is characterized by 
the presence of a large N V emission. However, compared to 
B 1246— 0542, the polarized spectrum is complex and difficult 
to interpret. In particular, the polarization does not increase in 
the Si IV trough while there are some absorption seen in the 
polarized flux close to the location of the Si IV emission line. 
This latter absorption possibly results from C IV absorption at 
velocities larger than 25000 km which could be superim- 
posed on the Si IV emission line. The Aim BAL is polarized 
as the continuum. The rotation of 9 near 4600 A is not signif- 
icant at the 2 a level, the uncertainty on 9 in this region being 
large due to a decrease of p. However, if real, it may also be 
due to complex absorptions that are difficult to identify in this 
rather noisy part of the polarized flux. 

3. B2225— 0534 This object is the prototype of the 

BAL QSOs with C IV P cygni type profiles. The spectropo- 
larimetry of B2225— 0534 has been published by Stockmann 
et al. ( IMlt . C ohen e t al. ( fT995l . Goodrich & Miller ( fT995l 
and Ogle et al. (1999). We confirm the main results found by 
these authors : the continuum is highly polarized with p ^ 5%, 
the broad emission lines are unpolarized with a large drop of 
p in the C IV emission line, and the polarization strongly in- 
creases in the broad absorption lines. The polarization angle 
is constant throughout the spectral domain. The emission lines 
are absent from the polarized flux except a small residual from 
Cm]. Since no indication of variability has been reported for 
this object by previous authors, the fact that our results are in 
very good agreement with theirs provides a useful cross-check 
of our reduction. 
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4. B2240-3702 This BAL QSO is quite similar to 

B2225— 0534. Indeed, it shows P Cygni type profiles and 
nearly black broad absorption lines. But B2240— 3702 is also 
different with very weak emission lines and a double trough 
structure in the C IV and Si IV absorptions. The polarization 
strongly increases in the low velocity part of the absorption 
reaching a maximum of ~ 25 % in C IV. There is a dip of the 
polarization at the position of the emission lines. The polariza- 
tion redward of C IV is flat (p ^ 2 %) with a small decrease at 
the location of C III] . The polarization angle rotates accross the 
broad absorption lines and is constant elsewhere. The broad ab- 
sorption lines appear in the polarized flux although shallower 
There is no evidence for any emission line in the polarized flux. 

5. B1331-0108 The data for this BAL QSO have a 

lower signal to noise ratio. There are however indications that 
the polarization increases in the C IV, Si IV and Al III absorp- 
tion lines while there is no clear evidence for a rise of p in the 
Mg II trough. Nothing can be said about the emission lines. 

6. B1413+1 143"^ This BAL QSO is a four component 

gravitational mirage (Magain et al. 1988 1. Again, the polariza- 
tion strongly increases in the absorption troughs. Although the 
absorptions of C IV and Si IV are nearly black in the total flux, 
they are shallower in the polarized flux. There is also a rota- 
tion of the polarization angle accross the C IV absorption line. 
Nevertheless, the most outstanding result in this object is the 
prominence of the emission lines in the polarized flux, indicat- 
ing that they are polarized at a similar level as the continuum. 
Only the semi-forbidden C III] line is absent in the polarized 
flux. Our results are in good agreement with those published 
previously by Goodrich & Miller (1995 1, Schmidt & Hines 
(fT999l and Ogle et al. (1999 1, despite of the fact that the to- 
tal flux and the polarization ofB1413+1143 are known to be 
variable (e.g. Goodrich & Miller fl995l Qstensen et al. ll997> . 

5.4. Definition of spectropolarimetric indices 

In this section we define four indices in order to describe the 
polarization of the emission and absorption lines in BAL QSOs. 
The measurement of these indices has been carried out on our 
spectropolarimetric data and on the spectra from Ogle et al. 
( fT999l . They are given in Table g] For the five BAL QSOs of 
our sample with good quality data and which were also ob- 
served by Ogle et al. (^999^ (i.e. B0059-2735, BlOl 1+0906, 
B1246-0542, B1413+1143 and B2225-0534), the measure- 
ments are done independently for both data sets. Since they are 
in good agreement within uncertainties, we only give average 
values in Table|3 

The first two quantities measure the polarization in the C IV 
and C III] emission lines relative to the polarization of the red- 
ward adjacent continuum. They are given as Pe/Pc CiV and 
Pe/Pc Cm] in Table 13 The other two indices are related to 
the polarization in the C IV absorption line. The first one is the 
maximal polarization in the C IV trough, Pmax- The second one, 
denoted SI, is the ratio of the C IV absorption depths (measured 
at the wavelength of the deepest C IV absorption) in the polar- 



Table 4. The spectropolarimetric indices. The measurements 
were done using the spectropolarimetric data illustrated in Figs. 
|3land0]and the polarization spectra of Ogle et al. il999t . The 
continuum polarization, pq, and its associated error, ap, are 
taken from Table 1. po, cTp and Pmax are in % while the other 
quantities are unitless. 



Object 


Po 


Op 


Pe/Pc 

Civ 


Pe/Pc 

Cm] 


SI 




B0019+0107 


0.98 


0.02 


0.90 


1.00 


- 


4.8+0.6 


B0043+0048 


0.09 


0.06 


- 


- 


- 


4.4+0.9 


B0059-2735 


1.49 


0.02 


- 


- 


- 


26.1+4.0 


B0105-2634 


2.41 


0.08 


0.50 


0.60 


0.56 


10.5+1.6 


B0137-0153 


1.09 


0.05 


0.50 


0.50 


- 


6.5+2.1 


B0146+0142 


1.24 


0.02 


0.40 


0.70 


0.60 


4.3+0.7 


B0226-1024 


1.81 


0.01 


0.60 


0.70 


0.50 


7.3+0.8 


B0842+3431 


0.51 


0.01 


0.80 


1.00 


0.60 


2.9+04 


B0903+1734 


0.67 


0.02 


0.30 


0.70 


- 


5.2+1.3 


T~i f\f\ ■~} ^ 1 C f\f\/' 

B0932+5006 


1.11 


0.02 


0.70 


0.80 


0.85 


2.7+0.7 


B 101 1+0906 


2.00 


0.08 


0.78 


0.66 


0.74 


9.0+0.8 


B1212+1445 


1.49 


0.03 




0.60 


1.00 


4.3+1.0 


B1232+1325 


3.19 


0.04 


0.50 


0.79 


0.88 


12.2+2.6 


B 1235+0857 


2.53 


0.07 




0.70 


1.00 


3.3+0.8 


B 1246-0542 


1.26 


0.01 


0.63 


0.77 


0.13 


7.7+0.6 


B1331-0108 


1.55 


0.14 


0.70 


0.30 




7.2+14 


B 1333+2840 


4.67 


0.02 


0.20 


0.50 


0.76 


8.5+0.9 


B1413+1143 


1.52 


0.04 


0.90 


1.00 


0.66 


13.2+2.2 


B 1524+5 147 


3.49 


0.01 


0.50 


0.70 


0.60 


6.0+0.2 


B2225-0534 


4.26 


0.02 


0.21 


0.70 


0.93 


10.7+2.9 


B2240-3702 


2.09 


0.19 


0.78 


0.74 


0.85 


25.5+2.1 



ized flux p X F\ and in the total flux F\ where the continuum 
has been normalized to unity {Fc is the flux in the continuum 
and pc is the polarization of the continuum) i.e., 



SI = 



P X Fa 

Pc X Fa 



1 - 



(1) 



Of course, we have < SI < 1. SI is a measure of the strength 
of the absorption in the polarized flux with respect to the ab- 
sorption in the direct flux. 

The values given in Table |4] are measured for 21 objects 
with good signal to noise data. The uncertainties on Pmax are 



Note that tiiere is an error in the polarization position angle data 
of B1413+1143 displayed in Lamy & Hutsemekers (il999j 



taken from Ogle et al. (I1999t or they have been measured from 
our data (Figs.|3land|3. The uncertainties on the other indices 
are more difficult to estimate, namely since they depend on the 
position of the continuum and the signal to noise in the po- 
larized spectrum. Based on two independent sets of measure- 
ments, we estimate that the typical uncertainties of these quan- 
tities are --10-15 %. 

5.5. Correlation analysis 

Table |5]gives the results of the correlation analysis performed 
between the four spectropolarimetric indices described in the 
previous section, the broad band polarization po and three im- 
portant and well-measured spectroscopic indices (BI, DI and 
as) given in Table 1. Correlations have again been searched 



Lamy & Hutsemekers: Polarization properties of broad absorption line QSOs : new statistical clues 



13 



Table 5. Results of correlation tests between spectropolarimet- 
ric indices, po and three indices derived from optical spectra. 
Pt and Pr, give the probabilities that the correlation between 
the two indices are due to chance, n is the number of points 
considered for the correlation. A correlation is emphasized if 
P < 0.05. The last column gives the sign of the correlation 



Index 1 


Index 2 


Pr 




n 




PO 


Pe/pc Civ 


0.067 


0.091 


17 


— 


Po 


P,/Pc Cm] 


0.147 


0.167 


19 


— 


Po 


SI 


0.224 


0.254 


15 


+ 


Po 


Pmax 


0.029 


0.019 


22 


+ 


Pe/Pc Civ 


Pe/pc Cm] 


0.012 


0.018 


17 


+ 


Pe/pc Civ 


SI 


0.440 


0.647 


13 


— 


Pe/Pc C IV 


Pmax 


0.966 


0.936 


17 


+ 


Pe/Pc C IV 


BI 


0.523 


0.583 


17 


— 


Pe/Pc C IV 


DI 


0.579 


0.640 


17 


+ 


Pe/Pc Civ 


OtB 


0.797 


0.763 


17 


+ 


Pe/pc Cm] 


SI 


0.553 


0.714 


15 




Pe/pc Cm] 


Pmax 


0.823 


0.777 


19 




Pe/pc Cm] 


BI 


0.970 


0.951 


19 


+ 


Pe/pc Cm] 


DI 


0.852 


0.913 


19 


+ 


Pe/pc Cm] 


CtB 


0.823 


0.775 


19 




SI 


Pmax 


0.959 


0.947 


15 




SI 


BI 


0.919 


0.914 


15 


+ 


SI 


DI 


0.033 


0.073 


15 




SI 


as 


0.203 


0.221 


15 


+ 


Pmax 


BI 


0.008 


0.004 


22 


+ 


Pmax 


DI 


0.080 


0.093 


22 




Pmax 


Ob 


0.089 


0.059 


22 


+ 



O ■ X 



2 4 6 8 



Fig. 5. The correlation between the maximum value of the po- 
larization in the C IV absorption line and the balnicity index 
BI. Open squares are HIBAL QSOs, filled squares are LIBAL 
QSOs and crosses are unclassified BAL QSOs. The open circle 
corresponds to an uncertain measurement (see text). 



for using the Kendall r and Spearman Ts rank correlation coef- 
ficients. With only 21 objects, we relax our criterion and con- 
sider as significant correlations with P < 0.05. With 22 corre- 
lations (assumed to be independent), we expect that one corre- 
lation with P < 0.05 is due to chance. 



0.4 0.6 
Pj/Pc CIV 



Fig. 6. The correlation between the polarization degrees in the 
C IV and C III] emission lines relative to the polarization in 
the adjacent redward continuum. The dashed line refers to an 
equal polarization in the two lines. Symbols are identical to 
those used in Fig.|5l 



Although the relation Pmax — BI (see Fig. |5} is the most 
significant correlation in Table|5] we must consider it with cau- 
tion. Indeed, Pmax being measured where the absorption is the 
deepest, it may depend on the spectral resolution, more partic- 
ularly when the residual intensity in the trough is nearly zero, 
which usually corresponds to large BI. In this case, it is possi- 
ble that the higher the resolution, the lower the flux in a single 
pixel and the larger the measured Pmax- However, a majority of 
BAL QSOs with large BI do not have completely black troughs, 
and the Pmax — BI correlation is most probably real. This cor- 
relation is discussed in more details in Sect. 6.2. 

As found in previous spectropolarimetric studies (e.g. Ogle 
et al. 19991, the polarization in the emission lines is systemat- 
ically lower than that the polarization in the continuum, which 
corresponds to pe/pc < 1 for both the CiV and Cm] emis- 
sion lines (see Table|4}. In addition, we find in Table|5]another 
significant correlation which relates the polarization in the C IV 
and C III] emission lines. As illustrated in Fig.|6l the polariza- 
tion in the C III] emission line is not only correlated to but also 
systematically higher than the polarization in the C IV emis- 
sion line. If the broad emission lines appear in the polarized 
flux, one possible explanation may be the effect of two dif- 
ferent polarizing mechanisms, one of them either depolarizing 



C IV or polarizing C III] . Lee (I1994t has shown that resonance 
scattering is more efficient at producing polarization in semi- 
forbidden lines such as C III] than in most permitted lines such 
as C IV. On the other hand, if there is only one polarizing mech- 
anism at work for both emission lines, the emitting regions may 
have different geometries and/or sizes with respect to the scat- 



tering region (e.g. Goodrich & Miller [1995). Cassinelli et al. 
fl987'( have calculated the polarization of a spherical emission 
region of radius R scattered by a region located at a distance 
D > R. They find that the polarization should be corrected by 
a factor y/1 — (R/D)^ compared to the case of a point source. 
If the size of the C III] emitting region is smaller than the size 
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Fig. 7. The correlation between the detachment index DI and 
SI, the ratio of the C IV absorption depths in the polarized flux 
and in the total flux. Both quantities are unitless. Symbols are 
identical to those used in Fig. [S] The open circle corresponds 
to an uncertain measurement (see text) 



of the C IV emitting region, the geometric dilution may induce 
a larger polarization in the C III] emission line than in the C IV 
emission line. Of course, one must keep in mind that neither the 
broad emission line region nor the scattering region are likely 
to be spherical. Since these results are in contradiction with 
reverberation mapping studies (Clavel et al. ^91 Chiang & 
Murray Ll 996 1, this explanation based on a purely geometrical 
dilution seems unlikely. In some BAL QSOs, the broad emis- 
sion lines do not appear in the polarized flux indicating that 
they are not polarized at all. In this case, the observed corre- 
lation may simply reflect the fact that the polarization in the 
continuum is larger in the blue part of the spectrum (i.e. at the 
wavelength of C IV) than in the red part (i.e. at the wavelength 
of Cm]). 

There is also a marginal {P ~ 0.03) correlation between 
the detachment index and SI, which is worth considering. In 
Fig.Q we plot SI versus DI. It clearly appears that if we do not 
consider the object represented by an open circle^, the corre- 
lation between SI and DI is much more significant. Indeed, in 
this case, the Kendall t and the Spearman statistics give a 
probability lower than 0.005 that these two quantities are cor- 
related by chance. Additional measurements are badly needed 
to confirm this possible correlation. The correlation between SI 
and DI indicates that in P Cygni type BAL QSOs, the BALs in 
the polarized flux are nearly as deep as in the direct flux (cf. 
B2225-0543), while in BAL QSOs with detached troughs, the 
BALs in the polarized flux are much shallower than in the di- 
rect flux (cf. B1246— 0542). As discussed in the next section, 
the consequences of such a relation for constraining BAL QSO 
models may be important. 



^ By carefully inspecting the data of Ogle et al. <1999S, we noticed 
that the C IV absorption line in B1212+1445 is partly cut at the lower 
end of the spectrum such that the value of SI given in Table |5lf or this 
object might be highly uncertain. 



6. Discussion 

6. 1. The principal component PC2 : tine 
anti-correlation po - DI 

The fact that po is only coiTelated with the detachment index 
emphasizes the importance of DI to understand the geometry 
of the outflow in BAL QSOs. The PCA has confirmed the im- 
portance of the pq — DI anti-correlation which is dominant in 
PC2, one of the two significant principal components. 

The WfD model of Murray et al. (|1995 1 proposes a natural 
explanation of DI, hence of the detached troughs observed in 
some BAL QSOs (see Fig.|8} . According to this model, an ob- 
server looking through the wind close to the disk (edge-on) sees 
absorption at all velocities including at low velocity. P Cygni 
type profiles are then produced and the observer measures a 
small value of DI. On the contrary, if the line-of-sight towards 
the observer grazes the upper edge of the wind, it misses the 
low velocity part of the absorption, the wind streamlines be- 
ing more vertical than radial at low velocities. Therefore the 
absorption may start several thousand kilometers per second 
blueward of the corresponding emission peak, resulting in de- 
tached troughs and large values of DI. In the WfD model, DI is 
naturally related to the orientation of the accretion disk with 
respect to the line of sight. On the other hand, in the torus 
model (Schmidt & Hines Il997t . DI has no simple explana- 
tion. In polar wind models such as the one described by Punsly 
(1999 1, the detached troughs are caused by intermittent varia- 
tions of matter ejections along the bipolar outflow. In this case, 
the line detachment is related to a time-dependent phenomenon 
and does not depend on orientation. 

Since the continuum polarization may also depend on the 
inclination, the relation between po and DI favours a BALR 
model which is orientation dependent. In the WfD model, the 
dependence of po with DI is naturally explained since higher 
polarization is expected close to the disk where attenuation of 
the direct unpolarized continuum is larger As the inclination 
(measured from the disk axis) decreases, the growing impor- 
tance of the direct continuum gives rise to a dilution of the 
polarization. Therefore, detached troughs are naturally associ- 
ated with lower optical polarization. Since the polarization is 
not correlated with the slope of the continuum, it rules out an 
extinction by small dust particles as being the dominant extinc- 
tion mechanism. A significant part of the attenuation should 
then be neutral, dominated either by electron scattering or by 
extinction by large dust particles. 

The scatterers at the origin of the polarization may lie in the 
inner parts of the wind itself or they may be located in a polar 
region, which is usually assumed to interpret the spectropo- 
larimetric observations (e.g. Hines & Wills 1 19951 Goodrich & 
Miller 1995 1, although its origin is unclear 

If the scatterers are located in the wind, the electron scat- 
tering model within the cyUndrical sector geometry depicted 
by Brown & Mc Lean ("1977'( may be considered to explain the 
Po — DI anti-correlation (see also Paper I). More particularly, as 
discussed in Sect. 3, Fig.^ shows that the relation po — DI is 
not a true anti-correlation : for low values of DI (P Cygni pro- 
files), there is a large range of po values. At a given inchnation, 



Lamy & Hutsemekers: Polarization properties of broad absorption line QSOs : new statistical clues 



15 



the Brown & Mc Lean J1977> model predicts a polarization 
which increases with the half-opening angle of the wind. The 
range in pq for a given inclination may therefore simply reflect 
a range of opening angles of the equatorial wind. For large val- 
ues of DI (profiles with detached troughs), the strong dilution 
by the direct unpolarized continuum leads to small values of po 
whatever the opening angle of the wind. The larger polarization 
of LIBAL QSOs compared to HIBAL QSOs (Paper I, Lamy 
|2003 1 may then be explained either by the stronger extinctions 
obser ved in LIBAL QSOs (Sprayberry & Foltz .1992 Reichard 
et al. I2003t or by a larger opening angle of their wind. This 
latter possibility is in agreement with the high covering factors 
needed to interpret their very weak [O III] emission (Boroson 
& Mevers FT^ Yuan & Wills 2003). This simple explanation 
qualitatively reproduces the observed pq — DI correlation. For 
HIBAL QSOs, the correlation is possibly weaker because of 
the variation from object to object such that the smaller range 
of Po values may mask the correlation. 

6.2. The anti-correlation SI - DI and the correlation 

Pmax BI 

If the scatterers are located in a polar region, part of the light 
from the continuum may be scattered and subsequently ab- 
sorbed in the equatorial wind. The variation of SI with DI, de- 
rived from the spectropolarimetric data, is naturally explained 
within this geometry, as illustrated in Fig.|8l: for nearly equa- 
torial line-of-sights (DI <C), the scattered flux crosses regions 
of the equatorial wind located close to the accretion disk where 
the opacity in the absorption lines is large. Therefore the po- 
larized flux is significantly absorbed and shows deep broad 
absorption lines (SI ^ 1, cf. B2225— 0534). For less inclined 
line-of-sights (DI ^), the scattered flux crosses regions of the 
equatorial wind where the opacity is much lower resulting in 
weaker absorptions in the polarized flux (SI ^). The quantity 
SI becomes ~ in the extreme case where the scattered flux 
is not absorbed at all (cf. B 1246-0542). It is more difficult to 
understand this relation if the scatterers are located in the equa- 
torial wind itself than in a polar region. The SI-DI relationship 
appears therefore important for constraining the geometry of 
the BALR, and is definitely worth confirming with additional 
spectropolarimetric observations. 

The other coiTelation derived from the spectropolarimetric 
data, Pmax — BI, may also be simply explained with such a 
model. If BI is large, the direct unpolarized flux is strongly ab- 
sorbed in the line resulting in a smaller dilution of the scat- 
tered component at these wavelengths. Therefore large values 
of Pmax are expected in the case of strong absorptions. 

6.3. A "two-component wind" model ? 

Although the po — DI and SI — DI anti-correlations can be in- 
terpreted assuming an equatorial wind issued from a disk, the 
existence of an additional polar region where scattering takes 
place is apparently needed, at least in some BAL QSOs. This 
additional region may simply consist of another part of the out- 
flow with a lower density and a larger amount of free electrons. 




Quasar 

" direct spectrum 

polarized spectrum 

Fig. 8. Schematic illustration of the SI-DI relationship. The 
equatorial wind, represented by the curved streamlines, 
emerges from a geometrically thin accretion disk located in the 
equatorial plane. A mirror of electrons is located in the polar 
region and scatters some radiation coming from the central re- 
gions along the line-of-sight. This polarized component is then 
subsequently absorbed in the equatorial wind. When the disk is 
seen edge-on, P Cygni type profiles are observed since the ob- 
server sees absorption at all velocities including the low ones 
near the disk. For less inclined line-of-sights, detached absorp- 
tion profiles are observed since the observer misses the low 
velocity part of the absorption. 



Since the properties of these possibly related regions are dif- 
ferent, we call hereafter this model the "two-component wind" 
(2CW) model, emphasizing that it is a natural extension or a 
variant of the WfD model. 

The 2CW model is supported by several theoretical models 
which have been developed in the last years (Proga et al. ll998l 
[19991 Proga l2000ll2003l Pereyra et al. 2004 1. These authors 
generalize the Murray et al. model, in particular by including 
a magnetic field distribution. The outflow may then be radia- 
tion or magnetically driven or both. The simulations show that 
different types of winds can be produced (from being nearly 
polar to nearly equatorial) by varying the magnetic field, the 
total luminosity or the ratio between the luminosity of the disk 
and that of the central compact object. In particular, one may 
qualitatively reproduce a 2CW model with a dense equatorial 
wind bounded on the polar side by a much less dense stream. 

Observationally, the rotation of the polarization angle 9 ac- 
cross broad absorption lines (a remarkable instance is given in 
B1413+1143, see Fig. 13 may provide evidence for the exis- 
tence of two scattering regions in some BAL QSOs. Indeed, the 
continuum may be polarized by scattering from a polar and an 
equatorial distribution of electrons (or dust) at the same time. 
The scattered flux from the polar region has a polarization an- 
gle perpendicular to the disk axis while the polarization pro- 
duced in the disk is parallel. The resulting polarization angle 
depends on the relative importance of both components. If the 
scattered flux from the polar region is re-absorbed in the equa- 
torial wind, its contribution is reduced at these wavelengths 
leading to the observed rotation of the polarization angle. On 
the other hand, there may be only one source of scatterers, and 
the rotation of 6 may result from another polarizing mecha- 
nism, for example resonance scattering since in P Cygni type 
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profiles the broad absorption is partly filled with resonantly 
scattered emission. 

Besides polarization, another indicator of geometry is the 
radio morphology. Recently, some observational evidence sup- 
porting the existence of much smaller and weaker radio jets in 
radio-quiet QSOs have recently emerged (Kuncic 11999V Since 
the axis of the accretion disk is usually thought to be aligned 
with the radio jets, the comparison of the position angles of 
radio axes (PA) and optical polarizations (6) can in princi- 
ple give clues on the location of the scattering region. If the 
broad absorption lines occur in an equatorial wind, BAL QSOs 
are quasars seen edge-on and we should detect a radio axis 
close to the plane of the sky. The optical polarization should 
be roughly parallel to the radio axis if the scatterers are located 
in the wind or roughly perpendicular to it if they lie in a po- 
lar region. Unfortunately, because of their weak radio fluxes 
and their large distance, there are very few BAL QSOs with 
a resolved radio structure. In the low redshift BAL QSO PG 
1700+5153, Kellermann et al. (.1994) measure PA - 145° 
for the radio structure observed at the VLA while Ogle et al. 
( 1199 91 and Schmidt & Hines (1999) measure a polarization 
position angle 9 - 53°. In the radio-loud BAL QSO LBQS 
1138-0126, Brotherton et al. (2002i observe a double-lobed 
radio structure with PA ^ 52° characteristic of an edge-on 
source in unified radio schemes. Their spectropolarimetry of 
LBQS 1138—0126 shows a continuous rotation of the contin- 
uum position angle from 180° in the blue to 150° in the red 
part of their optical spectra. In both objects, the optical polar- 
ization is then more or less perpendicular to the radio struc- 
ture, suggesting that the scatterers are located in a polar region. 
However, we must keep in mind that at the VLA resolution, the 
weak radio axis may be twisted by some interactions with the 
surrounding medium. Therefore, higher resolution images are 
needed to confirm the VLA results. Recently, Jiang & Wang 
( I2003> performed the first VLBl observations of three BAL 
QSOs. Among them, J 1312+2319 has a two-sided structure 
with a position angle PA ~ 59° indicating an edge-on orienta- 
tion. We have recently obtained the optical polarization of this 
object for which we find 6 ~ 166° (Sluse et al. 2004) indi- 
cating that the two position angles are roughly perpendicular. 
Although there are only very few radio observations of BAL 
QSOs, they apparently support a 2CW model, where scattering 
occurs in a polar region. 

In the 2CW model, the BALR is located in the equatorial 
plane. However, based on radio observations of FIRST BAL 
QSOs with the VLA, Becker et al. (2000) argue that BAL 
QSOs are not viewed from any particular orientation because 
they observed both steep and flat spectra within a sample of 
BAL QSOs discovered in the FIRST catalog. While, in BAL 
QSOs, it is not sure that the radio spectrum is a good indi- 
cator of the position angle of the radio axis as it is in radio- 
loud sources, the recent VLBI observations of Jiang & Wang 
(I2003> may confirm the claim of Becker et al. (2000 1. Indeed, 
while J 1312+2319 is seen edge-on, J 1556+3517 has a com- 
pact structure and a flat radio spectrum suggesting it is proba- 
bly seen pole-on. However, J 1556+35 17 has also very unusual 
spectropolarimetric properties (Brotherton et al. I1997> . possi- 
bly making it a peculiar object. 



6.4. The principal component PC1 : the correlation 
Bl -Few 

Until now, we have concentrated on PC2 which is most proba- 
bly driven by orientation. But what drives PCI in BAL QSOs ? 
Yuan & Wills (2003 1 have shown that z - 2 BAL QSOs fol- 
low the main correlations linking the properties of emission 
lines as depicted in the so-called Boroson & Green eigenvector 
(11992 ). In particular, BAL QSOs have the strongest Fe II and 
weakest [O III] emission among high redshift quasars. For non- 
BAL QSOs, these properties are usually interpreted in terms of 
accretion rates of gas around the central massive black hole 
(Boroson & Green 1992). When the accretion rate increases, 
the disk inflates because of the larger radiation pressure. The 
optical depth becomes higher and the optically thick gas pre- 
vents the UV ionizing photons to reach the narrow-line re- 
gion where [O III] forms. At the same time, the Fe II emis- 
sion strongly increases because of a larger amount of X-rays. 
In this scenario, BAL QSOs are quasars which accrete matter 
with a rate close to the Eddington rate (Yuan & Wills I2003t 
such that they display the most extreme Fe II and [O III] emis- 
sion properties. In the WfD model, it is possible that the larger 
the accretion rate, the higher the gas density in the wind and 
therefore the stronger the absorption. This could possibly ex- 
plain the positive correlation observed between BI and Fe II A 
2400, which dominates PCI. 

7. Conclusions 

With new broad-band polarization data and new measurements 
of indices on good quality spectra from the literature, we have 
performed a systematic search for correlations among a sam- 
ple of 139 BAL QSOs. We find six non-trivial correlations. We 
use a principal component analysis to isolate the main corre- 
lations in only two principal components : PC 1 , dominated by 
the BI — Fe II correlation, is possibly related to the accretion 
rate of matter around the central massive black hole, while PC2, 
dominated by the pq — DI anti-correlation, could be driven by 
orientation. 

We have also searched for correlations among indices mea- 
sured on spectropolarimetric data for a sample of 21 BAL 
QSOs. We find that the polarization in the C III] emission line is 
systematically larger than the polarization in the C IV emission 
line, and that the maximum of the polarization in the C IV ab- 
sorption line is higher for BAL QSOs with large values of BI. 
Another important result is a possible correlation between the 
detachment index DI and SI, a quantity which measures the 
depth of the C IV BAL in the polarized flux relative to that one 
in the total flux. These statistical studies outline the importance 
of polarization to disentangle the geometry of the outflows in 
BAL QSOs. We have also emphasized that BI and more par- 
ticularly DI are critical parameters to describe the BAL phe- 
nomenon. 

The anti-correlations po — DI and DI— SI are naturally ex- 
plained with a model where the continuum attenuation and the 
line absorption occur within a continuous equatorial wind ac- 
celerated from a disk. The direct continuum is also assumed to 
be scattered by a "mirror" of electrons located in polar regions. 
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possibly belonging to a polar wind. For edge-on line-of-sights, 
the direct unpolarized continuum is attenuated in the dense 
equatorial regions near the disk. Therefore, in BAL QSOs with 
P Cygni profiles, a large continuum polarization is observed be- 
cause the relative contribution of the scattered polarized com- 
ponent is large. The scattered flux is also absorbed in the equa- 
torial wind and a prominent absorption line is observed in the 
polarized spectrum. For a less incUned line-of-sight grazing 
the upper edge of the equatorial wind, the direct unpolarized 
continuum is less heavily absorbed and dilutes the continuum 
polarization. The scattered flux is also less absorbed in these 
lower opacity regions of the wind such that for BAL QSOs 
with detached absorptions, lower continuum polarization and 
shallow troughs in the polarized flux are observed. 

However, with polarization data only, we cannot unambigu- 
ously locate the scatterers. A convenient way to solve this prob- 
lem is to compare the polarization position angle with the ori- 
entation of the radio-axis. Unfortunately, only 3 BAL QSOs 
have been resolved at the radio wavelengths and have polari- 
metric data available. Interestingly, in these 3 objects, the opti- 
cal polarization angle is nearly perpendicular to the radio-axis, 
as expected from a model where scattering occurs in a polar 
region. 

Acknowledgements. H. Lamy thanks Prof. J. Lemaire and the Belgian 
Institute of Space Aeronomy for giving him the opportunity of per- 
forming observations in La Silla (Chile) in March 2002. The authors 
would also like to thank the referee for the careful reading. 

References 

Angonin M.-C, Vanderriest C, Remy M. et al., 1990, A&A 233, L5 
Arav N., 1997, in Accretion Phenomena and Related Outflows, 

lAU colloquium 163, ASP Conference Series, Vol 121, D.T. 

Wickramasingue et al. eds, 580 
Baldwin J.A., Stone R.RS., 1984, MNRAS 204, 347 
Becker R.H., White R.L., Helfand D.J., 1995, ApJ 450,559 
Becker R.H., Gregg M.D., Hook LM. et al., ApJ 479, 93 
Becker R.H., White R.L., Gregg M.D. et al., 2000, ApJS 538, 72 
Becker R.H., White R.L., Gregg M.D. et al, 2001, ApJ 135, 227 
Begelman M., de Kool M., Sikora M,. 1991, ApJ 382, 416 
Boroson T.A., Green R.F., 1992, ApJS 80, 109 
Boroson T.A., Meyers K.A., 1992, ApJ 397, 442 
Boroson T.A., 2002, ApJ 565, 78 

Brandt W.N., Comastri A., Gallagher S.C. et al., 1999, ApJ 525, 69 
Brotherton M.S., Tran H.D., van Breugel W. et al., 1998, ApJ 487, 
L113 

Brotherton M.S., van Breugel W, Smith R.J. et al., 1998, ApJ 505, L7 
Brotherton M.S., Croom S.M., De Breuck C. et al., 2002, AJ 124, 
2575 

Brown J.C., Mc Lean LS., 1977, A&A 57, 141 
CassineUi J.P, Nordsieck K.H. and Murison M.A., 1987, ApJ 317, 
290 

Chaffee F.H., Foltz C.B., Hewett PC. et al., 1991, AJ 102, 1627 
Chiang J., Murray N., 1996, ApJ 466, 704 
Clavel J., Reichert G.A., AUoin D. et al., 1991, ApJ 366, 64 
Clarke D., Stewart B.C., 1986, VA 29, 27 
Cohen M.H., ogle PM., Tran H.D. et al., 1995, ApJ 448, 77 
de Kool M., Begelman M., 1995, ApJ 455, 448 
di Serego Alighieri S., Cimatti A., Fosbury R.A.E., 1994, ApJ 431, 
123 



di Serego Alighieri S., 1998, in Instrumentation for Large Telescopes, 
Rodriguez Espinosa J.M. (ed), Cambridge University Press, 287 

Everett J.E., 2002, astroph/02 12421 

Foltz C.B., Chaffee RH. et al., 1989, AJ 98, 1959 

Francis P.J., Wills B.J., 1999, in Quasars and cosmology, ASP 
Conference Series 162, Ferland G. et al. (eds), 363 

Gallagher S.C, Brandt W.N. et al., 1999, ApJ 519, 549 

Glenn J., Schmidt G.D., Foltz C.B., 1994, ApJ 434, 47 

Goodrich R.W., Miller J.S., 1995, ApJ 448, 73 

Goodrich R.W., 1997, ApJ 474, 606 

Green PJ., Aldcroft T.L., Mathur S. et al., 2001, ApJ 558, 109 
Hall PB., Anderson S.R, Strauss M.A. et al., 2002, ApJS 141, 267 
Hartig G.R, Baldwin J.A., 1986, ApJ 302, 64 

Hazard C, McMahon R., Webb J.K., Morton D.C., 1987, ApJ 323, 

263 

Hewett PC, Foltz C.B., Chaffee FH. et al., 1991, AJ 101, 1121 
Hewett PC, Foltz C.B., 2003, AJ 125, 1784 
Hines D.C., Wills B.J., 1995, ApJ 448, 69 
Hutsemekers D., 1993, A&A 280, 435 

Hutsemekers D., Lamy H. & Remy M., 1998, A&A 340, 371 (Paper I) 
Hutsemekers D., Lamy H., 2000, A&A 358, 835 (Paper 11) 
Hutsemekers D., Lamy H., 2002, in Mass outflows in Active Galactic 

Nuclei : new perspectives, Crenshaw et al. (eds), ASP Conference 

Series 255, 207 

Hutsemekers D., Sluse D., Cabanac R. et al., 2004, in preparation for 

A&A 

Jiang D.R., Wang T.G., 2003, A&A 397, L13 
Kellermann K.I., Doeleman S., Wright M.C.H. et al., 1994, AJ 108, 
1163 

Korista K.T., Voit G.M., Morris S.L., Weymann R.J., 1993, ApJS 88, 

357 

Kuncic Z., 1999, PASP 111, 954 

Lamy H., Hutsemekers D., 1999, in Structure and kinematics of the 
Broad Line Region, ASP Conference Series 175, Gaskell M. (ed), 
125 

Lamy H., Hutsemekers D., 2000a, A&AS 142, 451 

Lamy H., Hutsem6kers D., 2000b, A&A 356, L9 

Lamy H., 2003, PhD thesis : Contribution a V etude de la polarisation 

des quasars de type BAL, University of Liege 
Lee H.W., 1994, MNRAS 268, 49 
Lee H.W., Blandford R.D., 1997, MNRAS 288, 19 
Low F.J., Cutri R.M., Kleinmann S.G. et al., 1989, ApJ 340, 1 
Magain P, Surdej J., Swings J.-P et al., 1988, Nature 334, 325 
Melnick J., Dekker H., D'odorico S., 1989, EFOSC, ESO Operating 

Manual 4, Version 2, ESO 
Menou K., Vanden Berk D.E., Ivezic Z. et al., 2001, ApJ 561, 645 
Moore R.L., Stockman H.S., 1984, ApJ 279, 465 
Murray N., Chiang J., Grossman S.A., Voit G.M., 1995, ApJ 451, 498 
Murtagh F, Heck A., 1987, Multivariate Data Analysis, D. Reidel 

Publishing Company 
Ogle P.M., 1997 in Mass Ejection from Active Galactic Nuclei, N. 

Arav et al. (eds), ASP Conference Series 128, 78 
Ogle PM., Cohen M.H., Miller J.S. et al., 1999, ApJS 125, 1 
0stensen R., Remy M., Lindblad PO. et al., 1997, A&AS 126, 393 
Pereyra N.A., Owocki S.R, Hillier D.J. et al., 2004, ApJ 608, 454 
Press W.H., Flannery B.P., Teukolsky S.A., Vetterling W.T., 1989, 

Numerical Recipes, Cambdridge University Press. 
Proga D., Stone J.M., Drew J.E., 1998, MNRAS 295, 595 
Proga D., Stone J.M., Drew J.E., 1999, MNRAS 310, 476 
Proga D., 2000, ApJ 543, 686 
Proga D., 2003, ApJ 585, 406 
Punsly B., 1999, ApJ 527, 624 

Reichard T.A., Richards G.T., Schneider D.P, Hall PB. et al., 2003, 
AJ125, 1711 



18 



Lamy & Hutsemekers: Polarization properties of broad absorption line QSOs : new statistical clues 



Schmidt G.D., Mines D.C., 1997, in Mass Ejection from Active 
Galactic Nuclei, N. Arav et al. (eds), ASP Conference Series 128, 
106 

Schmidt G.D., Mines D.C., 1999, ApJ 512, 125 
Scoville N., Norman C, 1995, ApJ 451, 510 
Serkowski K., Matthewson D.S., Ford V.L., 1975, ApJ 196, 261 
Simmons J.F.L., Stewart B.G., 1985, A&A 142, 100 
Sluse D., Mutsemekers D., Lamy M., Cabanac R., 2004, in preparation 
for A&A 

Sprayberry D., Foltz C.B., 1992, ApJ 390, 39 

Stockman M.S., Mier R.G., Angel J.R.R, 1981, ApJ 243, 404 

Stone R.PS., Baldwin J.A., 1983, MNRAS 204, 347 

Stoughton C, Lupton R.M., Bernardi M. et al., 2002, AJ 123, 485 

Surdej J., Mutsemekers D., 1987, A&A 177, 42 

Tolea A., Krolik J.M., Tsvetanov Z., 2002, ApJ 578, L31 

Tumshek D.A., 1988, in QSO absorption lines: probing the Universe, 

J.C. Blades et al. (eds), Cambridge, 17 
Veron-Cetty M.P., Veron P., 2000, A catalogue of quasars and active 

nuclei, 9th ed. Garching: European Southern Observatory (ESO), 

ESQ Scientific Report no. 19. 
Voit G.M., Weymann R.J., Korista K.T., 1993, ApJ 413, 95 
Wardle J.F.C., Kronberg PR, 1974, ApJ 194, 249 
Weymann R.J., Morris S.L., Foltz C.B., Mewett PC, 1991, ApJ 373, 

23 

Yamamoto T.M., Vansevicius V., 1999, PASJ 51, 405 
Yuan M.J., Wills B.J., 2003, ApJ 593, 11 



